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^ Abstract 
> 

^ ' We study in detail the application of renormalisation theory to models of cluster aggregation and fragmen- 
^ tation of relevance to nucleation and growth processes. In particular, we investigate the Becker-Doring 
equations, originally formulated to describe and analyse non-equilibrium phase transitions, but more 
recently generalised to describe a wide range of physico chemical problems. We consider here rate coeffi- 
cients which depend on the cluster size in a power-law fashion, but now perturbed by small amplitude 
random noise. Power-law rate coefficients arise naturally in the theory of surface-controlled nucleation 
and growth processes. The noisy perturbations on these rates reflect the effect of microscopic variations 
in such mean-field coefficients, thermal fluctuations and/or experimental uncertainties. In the present 
paper we generalise our earlier work that identified the nine classes into which all dynamical behaviour 
must fall (Wattis J A D W and Coveney P V 2001 J. Phys.A: Math. Gen., submitted) by investigating 
how random perturbations of the rate coefficients influence the steady-state and kinetic behaviour of the 
coarse-grained, renormalised system. We are hence able to confirm the existence of a set of up to nine 
universality classes for such Becker-Doring systems. 
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1 Introduction 



The purpose of the present paper is to demonstrate the robust nature of renormahsation methods in the 
theoretical description of nucleation and growth processes. We show that even in the presence of random 
perturbations the methods presented in our previous papers [[I^, [1^ provide correct asymptotic solutions. 



This paper thus extends and completes our earlier published work ITJ, |19|) wherein we derived a renor- 
mahsation procedure for the Becker-Doring equations and applied it to the case where the aggregation and 
fragmentation rate coefficients are of power-law form, given respectively by ar = ar^, br+i = br^, r being 
the aggregation number of a cluster (that is, the number of monomer particles within it). As previously 



described ITT], [I9|, these coefficients describe the rate at which a cluster of size r adsorbs monomer to grow to 
size r + 1, and the rate at which a cluster of size r + 1 sheds a monomer; this choice is well known to be ap- 
propriate for the description of surface- limited nucleation and growth processes [0, [1^. The Becker-Doring 



cluster equations are then 

Cf — J f — ]^ J f J J p — diip C\ bj<-^'\^ c^_i_ ^ ^ ^ 1 . X ^ 

where Cr{t) denotes the concentration of clusters of size r. Here we generalise our earlier work to consider 
rates given by the formulae 

ttr = ar^ + 5r, hr+i = br^ + Er+l- (1.2) 

in which the extra terms represent random perturbations to the deterministic rates, and which will be further 
quantified later on. While other authors have considered the effect of noise in such systems, the main thrust 
of such analyses has been to elucidate the temporal evolution of perturbations on cluster size distribution; 
see, for example, the work of van Dongen & Ernst [0, p^ . 

Renormahsation theory has been widely applied in the analysis of equilibrium phase transitions in statis- 
tical physics 0. In statistical mechanics, the basic idea underlying renormahsation theory is the transition 
from a microscopic to a macroscopic description of some phenomenon by the systematic filtering out of un- 
wanted degrees of freedom. In equilibrium phase transitions, near a critical point the system looks the same 
on all length scales and this physical insight can be translated into a set of transformations which leave the 
essential physical properties of the system unchanged - a procedure known as renormahsation. The widely 
used term "The Renormahsation Group" (RG) is technically inappropriate since the transformation loses 
information, and so is at most a semi-group, while the procedure comprises many different ideas and distinct 
methods, rather than being a formal monolithic edifice as the definite article would imply. Of more recent 
interest is the application of renormahsation ideas to non-equilibrium phenomena. While the physical mo- 
tivation behind the RG programme of coarse-graining microscopic models still seems appropriate to obtain 
the macroscopic properties, the complexity of far-from-equilibrium dynamics is such that in practice each 
specific system must be shown to be suitably scale-invariant. 

Of central interest in the present paper is the late-time asymptotic macroscopic behaviour of complex 
dynamical systems. Bricmont & Kupiainen have taken ideas from renormahsation theory together with 
asymptotic methods for the analysis of diffusive processes including nonlinear parabolic equations [|l], 0, 
1^ while Woodruff has recast multiple timescale problems using renormahsation ideas in ^2, 25, |23|, |24 



Woodrfuff's method allows the separation of equations for larger-scale phenomena and small-scale dynamics 
from a more general theory. Velazquez has recently used a renormahsation technique in an attempt to 
draw together the theories of Lifshitz-Slyozov coarsening and nucleation as modelled by the Becker-Doring 
equations |jl6| . 



We apply underlying concepts from renormahsation theory to study the Becker-Doring equations, which 
were originally formulated to study the kinetics of first-order phase transitions. They describe the stepwise 
growth and fragmentation of clusters in terms of the rates of the individual processes wherein monomer 
particles join or leave each cluster. The Becker-Doring equations have recently been subjected to more 
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conventional analysis using matched asymptotic expansions [T^. Rather than use the coarse-graining 
approximation, which emphasises the discrete nature of the equations, that analysis concentrates on the 



large-time limit where continuum approximations become valid . We have recently applied generalisations 



of these equations to a wide range of physicochemical processes, ranging from those involving surfactant self- 
assembly [P, through RNA-polymer-formation |0, |], |T8| to cement-setting p], 0]. In these studies coarse- 
graining procedures reduce large systems of equations down to lower- dimensional- "mesoscopic" -dynamical 
systems capable of theoretical analysis using standard techniques from the theory of differential equations. 
The coarse-graining contraction procedure summarised below is analogous to other renormalization methods 
used in statistical physics. 

There is some similarity between our methods and Woodruff's approach: we write the microscopic ag- 
gregation number r of a cluster as r = (A — l)n + 1 + k where n is of mesoscopic size and /c is a microscopic 
correction; we then aim to determine the problem on the mesoscale in a form which does not require us to 
simultaneously solve the microscopic problem. Thus microscopic detail is filtered out, but we are able to 
construct a simpler model which remains valid on larger scales. 



Preliminary results of this work were reported in an earlier publication [jrT|. In the present paper, 
the particular and physically relevant example of simple power-law rate coefficientss is generalised by the 
addition of small random perturbations. These perturbations influence the system's steady-states and large- 
time kinetics and their effects are studied here in detail. Following a brief recapitualtion of the model and the 
coarse-graining scheme underpinning our renormalisation procedure in the remainder of the present Section, 
the perturbations are introduced in Section ^. The central part of the present paper is concerned with an 
analysis of the effects which these perturbations have on our renormalisation procedure, and how such noise 
influences the contracted description of the model (Sections |]). The main issue at stake is the stability of 
the models under such minor perturbations of rate cofficients: our analysis leads directly to the identification 
of a set of nine generic classes of behaviour and true universality classes whose asymptotic behaviour is 
independent of all microscopic details (Sections |] & We suggest physicochemical scenarios that may 
correspond to the generic classes identified by our renormalsation analysis; future work, especially of an 
experimental nature, will be helpful in relating real-world nucleation and growth processes to these distinct 
universality classes. 



1.1 The Becker-Doring cluster equations 

In this section we give a basic outline of the Becker-Doring system of equations and their properties (con- 
sult [|TT], [T9| for more details). Let us start with a system in which a precursor chemical, P, spontaneously 
decays to form the monomer Ci, at some rate kf{p) where p = p{t) = [P] is the concentration of P. Further, 
we assume that this mechanism is reversible, with backard rate /cb(ci). The monomer is allowed to aggregate, 
with clusters growing and fragmenting according to the usual Becker-Doring cluster equations. Clusters are 
formed by two processes: either by the next smallest cluster size coalescing with a monomer, or by the next 
largest size losing a monomer. Only such monomer-cluster interactions are permitted in the Becker-Doring 
model of nucleation; cluster-cluster interactions are ignored. The system is thus governed by 

p = kh{p,ci)ci~ kf{p,ci)p, (1.1) 



oo 



ci = kf{p,ci)p- kb{p,ci)ci- Ji-J2Jr, (1-2) 



r=l 



— Jf—i — Jr^ ij" ^ 2), — afCfC\ — hj-j^\Cfj^\^ (1.3) 

where Cr{t) represents the concentration of clusters containing r monomers, and the constants a^, &r are 
aggregation and fragmentation rates respectively. 
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In this paper we consider the case for which the monomer concentration (ci) is held constant; thus the 
Becker-Doring equations we are concerned with are 

Cf = J-r—l — Jr^ iX — 2)) Jt ~ Oi^Cj.C\ — 6j._|_iCj._|_i , (1.4) 

with ci a given constant; we leave the more general formulation ( |l.l|) -( p]T3D for a future paper. The assumption 
of a constant monomer concentration corresponds to the case where the precursor chemical supplies monomer 
at a rate given by p = — Ji — X^J^i Jr-, so that ci = 0. This assumption is made in situations where the so-called 
"pool chemical approximation" is valid, namely where there is a large source of monomer species entering 
into the system at a rate which maintains the monomer concentration essentially fixed and independent of 
time. 

Given the rate coefficients a^, hr the partition function is defined by Qi = 1 and a^Qr = ^r+iQr+i- This 
generates the equilibrium solution Cr = QrCi, which is an equilibrium solution of the constant monomer model 
(|1.4|) , and also an equilibrium solution of the generalised system (|1.1|) - (|1.3| ) for the particular concentration 
of precursor chemical given by p = khCi/kf. In both cases the equilibrium solution corresponds to zero fiux, 
that is Jr = for all r, and the fiux from precursor to monomer is also zero {kfp — kbCi = 0). For certain 
choices of rate coefficients ar,br and certain monomer concentrations ci, the equilibrium solution will not 
decay to zero in the limit r ^ oo. In these cases, an alternative steady-state solution will be approached in 
the large-time limit. This solution is given by a constant nonzero fiux through the system, that is J,. = J 
independent of r. This condition yields the family of solutions 

which contains the equilibrium solution as the special case J = 0. The steady-state fiux J is determined by 
requiring the concentrations to decay to zero in the large-r limit, giving 

^ ^ /x] ^ ( ) 



1.2 The Becker-Doring system with power-law coefficients 

We now proceed to consider the renormalisation of Becker-Doring models in which the cluster rate coefficients 
are of power-law form, a dependence which is of immediate relevance to the description of surface-limited 
aggregation processes. We assume the rate coefficients for aggregation and fragmentation are respectively 

ar = ar^, br+i = br^, (1.7) 

so that the parameter 6 = aci/b is useful for classifying dynamical behaviour. The parameter p determines 
the variability of rate with cluster size, with p > implying that large cluster sizes have larger aggregation 
and fragmentation rates, and p < giving rates which decrease with increasing cluster size, the latter case 
is the less physically relevant, but is nevertheless also studied here for the sake of completeness. Typical 
values for p are p = 0, |, |, |, 1 for the examples of linear chain polymerisation, coagulation kinetics in 
two space dimensions, diffusion-limited coagulation in three dimensions, surface-limited coagulation in three 
dimensions, and branched chain polymerisation, respectively. Since a cluster's volume scales with aggregation 
number r, if we assume that clusters are spherical then their surface area scales with r^/^ and their diameter 
with r^/^, accounting for the presence of these exponents. More general exponents can be manifest in other 
situations 0. 
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The partition function Qj. is defined by Qr = {a/bY ^ as in the p = case considered earher E^. For 



< 1 the system approaches the equihbrium solution given by solving = 0, that is = 9'''~^Ci. Note that 
this solution is independent of p, although the way in which the equilibrium solution is approached depends 
on p. For 6 > 1 the equilibrium solution diverges at large r. Instead, for 6 > 1 the system approaches one of 
a family of time-independent solutions in which all fluxes are equal; Jr = J for all r implies 

Cr = 0'~^Ci {l-j'y , ] ,\ . (1.8) 

Since, for 9 > 1 the sum is convergent in the limit r ^ oo, the flux which gives the least singular behaviour 
in this limit is 

/oo 
^ k-PQ-^ . (1.9) 
k=l 

1.3 Coarse-graining procedure 

Following the general coarse-grained contraction with constant mesh size A in aggregation number (so that 
we only retain the aggregation numbers r = A„ = (ra — 1)A + 1), the kinetic equations reduce to 

in = Ln-l — Lni {t > 2), L„ = UnXnX^ — (1-10) 

a„ = TaA„aA„+i . . . aA„+i_i, = 7'6a„+i6a„+2 • • • &a„+i, (l-H) 

where the retained coarse-grained cluster concentrations are relabelled with xi := Ci the monomer 

concentration; ci is not involved in the coarse-graining since it has a special role in the Becker-Doring 
theory |P, |ll|. The parameters are the coarse-grained aggregation and fragmentation rates, now 



representing the addition or removal of A monomers to or from a cluster (rather than just a single monomer 
which occurs in the full Becker-Doring system). This flux of matter is denoted by L„. The concentration 
Xn{t) is representative of the concentrations for cluster sizes (A„_i -f 1) < r < A„. The factor T represents 
a change of timescale which ensures that the large time asymptotic behaviour of the reduced system coincides 
exactly with the original fine-grained system in the case of size-independent aggregation and fragmentation 
rates (a^. = a, K = b). Technically speaking, the need to redefine the timescale makes this procedure a 
dynamical renormalisation. 

1.4 Coarse-graining of power law coefficients 



If the rate coefficients in the original formulation in eqns (|1.4| ) are determined by simple power laws, namely 
ar = ar^ and 6^+1 = br^, then the coefficients in the reduced model are 



Thus 



= {[(r - 1)A + l][{n - 1)A + 2] . . . [nA]f . (1.12) 

loga„ = A loga + log(nA — A + j) ~ X\oga + p \og{nX — X + x)dx 

= Aloga+pA [log(nA)-l + (l-n)log(l-i)] . (1.13) 

For large n this asymptotes to loga„ ~ A log a + pAlog(An), so for simplicity we shall take a„ = {aX^rfY, 
which differs slightly at small values of n. The backward rate coefficient is then /5„+i = [bX^n^Y. 
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The new system has its own 6 parameter determining the balance between aggregation and fragmentation 
rates in the system which, for the moment, we shall call 9 = arX^/(3r+i = 9'^, so contraction of the system 
maps 6 to 6^. The parameter 6 thus plays an important role in the renormalisation procedure, the fixed 
points of this mapping corresponding to 6 = 0, 1, cxo; hence such systems are of special interest to us. Also 
the contraction maps coefficients with exponent p to those with exponent pX. Thus following a contraction, 
there are only three limits to consider: small p (namely p = 0) and large p (positive and negative). 

The effect of coarse-graining a Becker-Doring system is to modify the rate coefficients, by the map 
{6,p) \—>- {6^,Xp). If A is allowed to take on large values, there are only nine combinations of {6,p) which 
merit attention, namely all possible combinations of 6 = {0, 1, oo} and p = 0,p > 0,p < 0. These nine cases 
and their associated fixed points will be the basis of the analysis presented in Section 3 of this paper, which 
follows on from Section ^ where we develop the general theory for perturbed rates for the coarse-grained 
Becker-Doring system with arbitrary A. 

2 General analysis of the role of noise in coarse-grained Becker- 
Doring systems 

2.1 Form of rate perturbations 

Whilst a simple power-law description of rate coefficients might give the correct general behaviour over a 
large range of aggregation numbers, in any given physical system there will be some minor deviations from 
such a regular descriptions which may be due to a number of separate effects. For instance, the presence of 
microscopic fluctuations and/or experimental uncertainties suggest that one should be concerned about the 
stability of such theoretical models to minor modiflcations in the assumed rate coefficients. Indeed, Figure 2.3 
of Lewis's article in the monograph |T5[ shows that the dependence of free energy on cluster size can in reality 



be much more complicated than the simple smooth curves obtained from simple derivations (for example by 
considering surface and bulk energies and assuming all clusters are spheres, discs or needles; or by working 
back from simple laws for aggregation and fragmentation rates). In mathematical terms, therefore, we wish 
to know whether small random perturbations to the rates alter the results we have derived previously using 
renormalisation methods. The hope is that they will not have a major effect, so that the renormalisation 
scheme is stable with respect to such perturbations and points to the existence of true universality classes. 
To analyse the effect of such perturbations, we perturb the rates from a simple power law form by small 
random amounts, in such a way that all rates remain positive. Thus we assume that the rate coefficients 
have the form 

ar = ar^ + 6r, K+i = br^ + Er+i, (2.1) 

where Sr,Sr are independent random variables, with small mean and small variance, so we can assume that 
6r ^ Er ^ 1 for all r. Many of the approximations used in the ensuing analysis rely on the noise being 
small {u <^ 1). Even with such a restrictive assumption, by calculating the first higher order term in the 
asymptotic expansion, many important and interesting results are obtained. Later on, we shall discuss cases 
in which u ^ 1 and u = 0{1). In the nine special cases analysed below, the constants a, b in eqns ( p^.l| ) will be 
taken to be either zero or unity. When limr_»oo ctr or lim^-^oo br is zero, the corresponding perturbation {6r, Sr 
respectively) will be assumed to be strictly positive and have mean v. Specifically, we shall assume that they 
are distributed according to some continuous probability distribution functions fs^x), /e(x) respectively, with 
/(■) = whenever 6,e lies more than O(z^) away from zero. The expectation operator E[-] is defined by 

E[g{5)] = J g{x)fs{x)dx, and E[g{e)] = J g{x)f,{x)dx, (2.2) 
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the integrals being over all real x. The variance operator V[-] is defined by 

W[g] = E[g'] - E[g]'. (2.3) 

The conditions on the mean reduces to E[5r] = and/or E,[6r] = v for all r; in cases where lim^^oo o^r and/or 
lim^^oo &r are unity or larger 5^, will be allowed to take negative values as well as positive and we assume 
that they have zero mean and their variance is 0[y^\ (e.g. E[5,.] = with ]E[5^] = (9(z/^)). In all cases the 
perturbations are assumed to be small enough that the total rate constants in equation ( p.l|) are all positive. 

2.2 Effect of noise on cluster partition functions 

We first derive and analyse the modified equilibrium and steady states which are approached in the large-time 
limit. The kinetics by which these states are achieved is described in detail later (Section |2.3| ). First we aim 
to find the effect which the perturbations to the rate coefficients (5^, ^r) have on the cluster partition function 
Qr in each of the nine special cases described above. Since the partition function is defined in terms of the 
rate constants, the presence of non-zero perturbations briSr will affect Qr- Due to the definition c\ = 1, the 
partition function is identical to the equilibrium configuration. However, in Cases III, VI, IX and V (if p > 1) 
it is not the equilibrium solution which is approached in the large-time limit; instead the system evolves 
to a steady-state solution. Perturbations to the rate coefficients (5r,£r) thus modify the steady-state and 
equilibrium solutions through the partition function, and such modifications will be found here. In general 
the partition function is given by 

Qr = n 7^, (2.4) 

k=l 

whilst steady-state solutions are determined by the constant flux condition J = Jr = flrCr — K+iCr+i, for 
all values of r. We shall pay particular attention to the large aggregation number (r) behaviour, since this 
determines the observed behaviour of the system on the mesoscopic scale, and is the region we shall probe 
later with the coarse-grained contraction of this perturbed Becker-Doring system. 

2.3 Effect of noise on cluster growth kinetics 

The equations we study in this section are the kinetic equations with noisy coefficients, namely 

= (a(r - ly + Sr-i)cr-i - (arP + 5^)c^ - (6(r - 1)^ + e^)c^ + (for^ + e^+i)c^+i, (r > 2). (2.5) 

Having found the states which are approached in the large-time limit, we now find time-dependent solutions 
by using large-time asymptotics. We write the solution as Cr{t) = Qrip{r,t) or Cr(t) = cf^ilj{r,t) in order to 
investigate the manner in which il){r,t) 1 ast ^ +oo. In general, this leads to 

Ipr = ar{lpr+l " A) " KiA " ^r-l) " " ^r-l), (2.6) 



where J represents the steady-state fiux into the system; convergence to equilibrium corresponds to J = 0. 

2.4 Effect of perturbations on rate coefficients in the contracted system 

The rate coefficients following contraction are given by 

A„+i — 1 A„-|.i — 1 

an = T Yl «r-, Pn+1 =T J| br+l- (2.7) 

r=A„ r=A„ 
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where we shall use T = A so that the 'clean' coefficients = ar^, = br^ are mapped to a„ = an^'^, 
Pn+i = PnP^ with a = a'^ and (3 = (see Section |1.4| and [0). Since the nine special cases which will be 



studied in detail later correspond to a, h equal to zero or unity, following the coarse-graining contraction we 
are concerned with establishing like-for-like correspondence in which a, (3 equal zero or unity. Here T is chosen 
to simplify the algebra; an alternative expression could be used, with a consequent increase in the complexity 
of the ensuing equations. Alternative choices for T merely influence the units of time, and have no effect on 
the equilibrium or steady-state solutions or the effect of the rate perturbations; the large-time asymptotics 
will only be affected by a linear transformation to the time variable, t. The problems we are concerned with 
here involve properties of the noise following contraction, for example whether its amplitude is dependent on 
aggregation number, and its order of magnitude as a function of v. We now reduce the randomly perturbed 
coefficients from ( p.l|) in the microscopic description of cluster-formation to the mesoscopic description in 
which an = avP^ + A„, [3n+i = PnP^ + En+i with the aim of finding how the perturbations An, En in the 
contracted model depend on the perturbations 6r,er in the full description of the model. In some cases this 
is trivial, since a = or 6 = 0, and in these cases A„, En+i is the product of a set of Sr,£r, but in other cases 
the relationship is more complex, and in such cases we give the leading-order expression for A„, En+i in the 
asymptotic limit u —>■ 0. We have 




On = J I n "'''"^ 




1+^ 

arP 



'A„+i-l \ / A„+i-l c 

^1 n H 



since 6r = 0{v) <^ ar^ = 0{1)] correction terms to the above approximation are thus O^u"^). Using ( |1.12| )- 



(|1.13 ), the leading-order term is approximated by avP'^ where a = Ta^X^^ and the ffist correction term due 



to the rate perturbations is defined by 



An+1-1 r 

An = anP^ E (2.9) 



r=An 



Similarly for the fragmentation rates we have Pn+i = l^fT-^^ + -E'n+i with (3 = Th^X^^ and 



An+1-1 



r=A„ 



En^,=f3nP^ E (2.10) 



2.5 Consistency and accuracy analysis 

Having derived coarse-grained rate coefficients which include the leading order corrections due to random 
perturbations from a power law, we now use these rate coefficients to construct a partition function T„ for 
the coarse-grained model from the macroscopic rates This is determined by 



which implies 



Ti = 1, a„T„ = /3„+iT„+i, {n > 1), (2.11) 

an ^'^ 
'^Af = n ^"^^ logTjv = ^ (loga„ - log/3„+i). (2.12) 

n=l Pri+1 n=l 

We can then compare this result with the full partition function as a check on the coarse-graining procedures, 
with the hope of finding T^v = Qr when i? = A^v to leading order in v. For cases in which a steady-state 
solution is approached as t ^ 00, we compare the steady-state solution x^^^ calculated from the noisy 
contracted rates, with the steady-state solution of the full model, cf^. If moreover Cr = Xn when 

r = A„ then we will have demonstrated that no crucial information has been lost in the coarse-graining 
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reduction. This test thus amounts to a consistency check which we shall carry out not just at leading order, 
but also to first order in v (that is in A^jE'^+i) for the purposes of accuracy assessment. These checks will 
enable us to demonstrate if and when we have successfully bridged the scales from microscopic to mesoscopic 
for certain choices of the coarse-graining parameter A. 



2.6 Effect of noise on the late-time kinetics of the contracted system 



Having checked that the equilibrium and steady-state solutions of the microscopic description of the system 
have been faithfully reproduced in the coarse-grained system, we perform one final calculation to check that 
the late-time asymptotic form of the kinetics is also correctly replicated. In general the contracted system of 
equations can be written using ( p..lO| ) and (p78| ) - (|2.10|) as 



a\n- 



*^n— 1 



/3(n-l)P^ + E, 



an 



+ A„ 



PnP^ + E, 



n+l 



(2.13) 



In the following section we examine the equilibrium and steady-state solutions of this model, as well as 
the manner in which the time-dependent solutions convergence to these solutions. We then compare this 
behaviour with the properties of the microscopic system prior to coarse-graining. We now analyse in greater 
detail the nine cases summarised at the end of section [L^; that is those corresponding to p = 0, p > and 
p < in the three cases ^ = 0, 1, oo. The conditions 6^ = 0, 1, oo correspond to (a,/3) equalling (0, 1), (1, 1), 
(1,0) respectively. The fourth case (a,/3) = (0,0) yields a special limitting system which will be described 
in Section |3.10| . 



3 The nine generic classes of asymptotic behaviour 

Nine distinct generic classes of asymptotic behaviour were picked out by the RG analysis in Section 1.4 (see 
also []TT], ^).For each class or case we shall perform the five calculations described above. We start with the 
effect of the perturbations on the partition function, and hence the equilibrium solution. In cases where a 
steady-state is approached rather than the equilibrium we also calculate the modified form of the steady-state 
solution. We then calculate the effect on the kinetics of approach to steady-state or equilibrium. Following 
these two calculations we turn to the contracted description of the system. The latter three calculations 
combine the problems of dealing with the noise and the coarse-grained system. We analyse the contracted 
Becker-Doring system derived from the full microscopic model with rate coefficients given by a combination 
of a power-law with random perturbations, and investigate whether coarse-graining averages noise out of the 
system or exaggerates its effect. In so doing, we discover how perturbations to the microscopic rates manifest 
themselves in the rates of the contracted system. We then establish consistency by showing that the partition 
function calculated from the coarse-grained rates is the same, to leading order, as the partition function from 
the microscopic system sampled over the coarse mesh. Finally, we analyse the effect of the perturbations in 
the mesoscopic model on the large-time kinetics of the coarse-grained system. 

3.1 Case I 

This case corresponds to a fragmentation-dominated system where the fragmentation rate is independent of 
cluster size, that is the rate at which large clusters shed monomers is the same as that for small clusters. 
Such situations may arise in certain types of polymer degradation, if the breakdown of chains takes place 
through end-monomer fission processes. 
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In the non-perturbed version of this case, = for all r, so the partition function is identically zero. 
Introducing noise to the aggregation coefficient makes the partition function non-zero 



r — 1 r /) — 1 \ / r— 1 



k = l 



a = nT^^-(n4) (i-E^^+i+^(^'))- (3-1) 

1 + Ek+l Vfc=l / V k=l I 



Thus the partition function depends primarily on the perturbations 5^; the perturbations Sr only enter at 
higher order. Since ¥\5r\ = u <ti 1, Qr decays rapidly to zero with increasing r, according to Qr = Oiy"^'^), 
thus Qr+i -C Qr for all r. Case I may hence correspond to a situation in which finite-length polymers or 
oligomers break down via stepwise fission processes. 

In the noisy case, there is a non-zero partition function Qr to which the system converges. To determine the 
manner of this convergence we transform from Cr{t) to new coordinates, Cr = Qr'4^r{t) so that at equilibrium 
ipr = '^- Then V'r(^) satisfies 

Tpr = A-l - V'r + Srilpr+l " Ipr) " ^rii'r " i'r-l)] (3.2) 

for large aggregation numbers (r) and at large times (t), the continuum limit is valid and this equation goes 
over into the following partial differential equation: 

^4||a,,.,,„)_^a,,_,,. ,3.3, 

The solutions of this type of equation typically approach equilibrium via a diffusive wave, whose position 
we shall denote by r = s{t). We transform from r to the new position variable z = r — s{t) relative to the 
wavefront. The leading order terms are then those involving formally yielding the equation s = l + Eg — Ss 
for the position of the wave. The expected value of the first perturbation (due to 5r) is zero, since E,[er] = 
for all r, but that of the second (due to Er) is positive, that is K[6r] = v for all r. Thus we have E[s] = 1 — z/, 
and the presence of perturbations slows the wave. This is not, however, a leading order effect since small noise 
makes only a small difference to the system's approach to equilibrium. Higher order terms from equation 
(|3.3|) lead to a description of the shape of the wavefront. If we start the system from a state of compact 
support (that is Cr(0) = for r > R for some R < oo), then the large-time and large aggregation number 
asymptotic solution is given by 

c,(t) = iQ,erfc|^^^ii^l, (3.4) 
'2{l + u)t^ 



where Qr is the partition function defined by equation (|3.1| ). Thus as well as travelling slightly more slowly 
than in the noiseless case, the wave is widened slightly, being spread out over a larger number of cluster sizes. 
These effects would be hard to measure experimentally since the typical concentration of large cluster sizes 
would be extremely small. Both these effects are due to the differing expected values of the perturbations to 
the forward and backward rate constants (namely that E[er] = whilst E,[6r] = v). 

Moving now to the coarse-grained, mesoscopic description of the problem, we find = 0, 6^ = 1 iniphes 
an = An and Pn+i = 1 + where 

A,>fi— 1 A^i— 1 

A„= n ^r, En+1= E Sr+1 + 0{U^). (3.5) 

r=A„ r=A„ 

Thus following an application of the coarse-graining contraction procedure, to 0{1) the rates are a„ = and 
Pri+i = 1 with small perturbations to the rates, = 0(i^^) and Er+i = Oiv). Thus the perturbations to the 
aggregation rates are much smaller than those to the fragmentation rates, making the system appear more 
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fragmentation-dominated than the uncontracted version. However, this effect is due to the contracted model 
only representing some cluster sizes: growth from Xr to Xr+i requires the addition of A monomers and hence 
happens at a rate of (9(z/'^) rather than the 0{i>) rate for growth from to c^+i in the full model. Thus the 
contracted description is dominated by fragmentation to the correct extent. Since the perturbations in the 
full model are independent, and distributed with E[(5r] = i^, E[£r] = 0, Y[5r\ = 0{i''^), Y[er] = 0{u'^), in the 
contracted description we have E[A„] = v^, E[E„] = 0, V[A„] = 0(1/^^), and N[En] = 0{iy^). 

We now use the contracted rates given by a„ = A„ and Pn+i = 1 + En+i to calculate a partition function 
for the coarse-grained system. Using ( |2.12| ) and ( ^.5| ), we find 

/N-l \ / N-1 \ /Ajv-1 \ / Aiv-1 \ 

= n 1 - E ^"+1 + (^i^') = n 1 - E ^-+1 + ^(^') • (3-6) 

\n=l / \ n=l J \ r=l / \ r=l / 

Thus the first two terms of the equilibrium solution of the reduced model agree exactly with that of Qn with 
i? = in the full model Q. 

In this case we know that the system will tend to the equilibrium solution given by x„ = T„ as in ( p.6|) , 
which closely approximates Qr ( |3.6| ) with r = A„. To find the way in which the solution is approached, 
we use a similar method to that successfully applied to the miroscopic description above, namely that of 
transforming to new dependent variables ipn{t) given by a;„(t) = Tnipn{t)- Then to first order in u we have 

1pn = 1pn-l-i'n- Enilpn-1pn-l), (3.7) 

since in this case A„ = Olu^) and = Olu). For large times and large ra- values the continuum approxima- 
tion 

^..,,.„)0_,,.„)|J. ,3.8) 

is formally valid. Since E[ii^„] = 0, the diffusive wave travels at unit expected speed and suffers no O(z^) 
correction term. Thus if the system is initiated from a;„(0) = for n > 2, the equilibrium is reached is via 
a diffusive wave moving from n = 1 to large n leaving the new equilibrium solution behind it as described 
by Xn{t) ~ |T„erfc((n-t)/v/2t). The noise thus has no effect at leading order, or at 0{i>), although higher 
order terms will influence the evolution of the system. Thus, in this case, coarse-graining has reduced the 
effect of the noise, since in the non-coarse-grained case, the kinetics are affected by Olu) terms, decelerating 
and widening the wavefront as described by equation 



3.2 Case II 

In this case the aggregation and fragmentation rates are both relevant and both size independent. This 
situation can be expected to arise in a wide variety of nucleation and growth problems, for example it is 
likely to pertain for cluster formation at saturation or low supersaturation levels during crystal growth. 

In this case introducing perturbations Sr, Br to the rates = ar^, 6^+1 = br^ as in (|2.1j ) modifies the 
partition function from Qr = 1 for all r to 

~ 1 + - Sk+i) + 0{u'). (3.9) 

k=l 

Small amplitude noise in the coefficients thus does not affect the leading-order behaviour of the system, at 
small values of r. To examine the large-r behaviour, we use the central limit theorem. For simplicity we 
assume that each of the random variables Sk,ek is distributed uniformly on the interval [—u,!/], thus each 
has an expected value of zero and variance of = ^v"^ . The difference 5^ — Sk+i thus has mean of zero 
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and variance of |z/^; and the above sum (equivalent to Qr — 1) has mean zero and variance \{r — Vju"^. At 
large values of r, the central limit theorem implies the sum can be approximated by a normally distributed 
random variable with zero mean and variance cr^ = |rz/^. Thus perturbations have a cumulative effect and 
may become significant when r = 0{v~'^)\ at this order of magnitude, the approximation in equation ( |3.9|) 
ceases to be valid, since neglected higher order terms then become significant. 
For Case II the kinetic equations are 

= (1 + 5,._i)Cr_l - (1 + £r)Cr - (1 + ^r)Cr + (1 + ^r+l)Cr+l, (3.10) 

where 5ri£r can be positive of negative, and are 0{v) with E[(5r] = IE[£:r] = 0. Formally, for large-times and 
large aggregation numbers, this leads to the continuum equation 

f-(l + K*.+^.))|| + (^..-.)f. (3.11) 

for ip{r,t) = Cr(t)/Qr. Since W.[6r — Er] = 0, and E[((5,. — e^)^] ~ <^ 1, diffusion dominates the advection 
terms, and there is no overall driving force on the diffusive wave. Effectively, it is pinned at r = 1 and 
equilibrium is reached by purely diffusive mechanisms. If we use compact initial conditions (that is Cr(0) = 
for all r > i? for some R < oo) then the large-time asymptotic solution is 

Cr{t) = Q.erfc (^^^ . (3.12) 

Note that E[(5r] = = K[6r] implies that the width scale is unchanged by the perturbations, in contrast with 
Case I, where the wave was widened due to E[5r] > 0. 

Following the coarse-graining contraction, the reaction rates are given by = 1-|-A„ and Pn+i = l+-E'n+i, 
where 

A„4-l— 1 An+l— 1 

^n= 5r + Oiu'), En+i= E er+i + Oiu^). (3.13) 

r=A„ r=An 

Thus a system in which aggregation and fragmentation are balanced is mapped to a similar system following 
a coarse-grained reduction. Since K[6r] = = E[£:.r] and ¥[5,.], ¥[5^] = O^u"^) we have E[A„] = = E[i?„] and 
V[A„], V[£'n] = Olu). The central limit theorem implies that the variances of the noise will increase linearly 
with A. Since both 6k, Ek have variance proportional to z/^, Ar,Er have variance proportional to Az/^. So in 
order for the noise to remain a small correction term, one has to ensure that A ^ z/^^; for small perturbations 
(z/ -C 1) this does not constitute a significant restriction. 

Using ( p.l2| ) and (|3.13|) , we find the partition function T„ for the coarse-grained system is given by 

Af-l A]v-1 

= 1 + 5: (A, - + 0(z.2) = 1+ i^r - Sr+l) + 0{u^). (3.14) 

n=l r=l 

Thus the first two terms of the equilibrium solution of the reduced model agree exactly with that of Qr with 
_R = A^r in the full model (|3.9|) . Thus in the presence of noise, the contracted system tends to an equilibrium 
solution of the same form, to first order, as the full system. 

The kinetic equations determining the approach to equilibrium are 

Xn = {l + An-l)Xn-l - (1 + En)Xn " (1 + A„)x„ + (1 + K+O^n+i, (3.15) 

where An, En are 0(z/) with E[A„] = E[i?„] = 0. The system tends to the equilibrium solution x„ = T„ by 
a purely diffusive mechanism as described formally by the continuum limit equation 
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for ip{n,t) = Xn{t)/Tn, so that ip ^ 1 as t —>■ oo. Since the expected values of perturbations A„ and 
En are both zero, there are no Olu) correction terms to ip{n,t) in the large-time limit, and we have the 
solution Xn{t) = T„erfc(?2/2-\/t). Thus to O^u) the large-time kinetics of the system are almost identical 
to the uncontracted model. Using r = A„, Cr = x„ and Qr = T„, the above formula for x„(t) implies 
Cr = Qrerfc(r/2Av^) whereas analysis of the full model gave (|3.12| ); thus the only difference between full and 
contracted models is in the timescale. 



3.3 Case III 

This case models situations which favour the formation of clusters since the system is dominated by aggrega- 
tion; however the rate of aggregation is size-independent so that nucleation and cluster growth processes are 
balanced. One physicochemical scenario that this Case may describe arise in certain surfactant self-assembly 
processes, in which amphiphile monomers attach to a growing assembly (for example, wormlike micelles, 
vesicles, and so on). 

In the pure aggregation case with no noise the partition function is not defined since all the fragmentation 
coefficients bj. vanish. The presence of noise makes it possible to define the partition function 

Qr = Cli f 1 + E 4 + Oiu')] . (3.17) 

This is strongly dependent on the perturbations e^, implying rapid growth in with r, specifically Qr = 
0{i'^^'^^^^). However, the system does not approach the equilibrium state Cr = Qr, rather it tends to a 
steady-state configuration. In the unperturbed problem this state is = 1 for all r, which has the steady- 
state fiux J = 1; when noise is added to the rate coefficients this state is modified to J = 1 + (5i— £2) + C>{i''^) 
by the calculation (|1.6|). Since the expected values of Sk,ek satisfy E[5fc] = and ¥.[ek] = z/ > 0, the expected 
value of the steady-state fiux is reduced by the presence of small amplitude noise from J = 1 in the noise-free 
system to E[J] = 1 — u. The concentrations asymptote to the modified steady-state 

C = 1 + (5l-£2 + £r-+l-5r) + Oiu^). (3.18) 

The presence of Olu) noise in the reaction rates alters the steady-state solution by an amount of Olu). Note 
that the noise in the two rate constants ai, 62 affects the limiting concentrations of clusters of all sizes. 

The kinetics of this case are governed by the approach to the steady-state solution (not the equilibrium 
solution). Since the perturbations to the fragmentation rates (e^) are all positive and the perturbations to the 
coagulation rates can be positive or negative, we have ^[Sk] = and ^[ek] = We replace the differential- 
difference system by a partial differential equation by taking the continuum limit and then seeking a diffusive 
wave solution of this partial differential equation. First, we transform variables by Cr{t) = cfj^^^iprit) so that 
'?/'r(t) ^ 1 as t — >• 00 according to 

lpr= (1 + 5r)(^r+l " A) + Srilpr-l " A) + ^(^r-1 " ^r+l), (3-19) 



«sss 



which formally goes over to the continuum limit 

^ = 1{1 + Sr+ Br)^ -{l + Sr+Sr- 2e.+l)^. (3.20) 

We define s{t) to be the position of the wave and transform from r to an as yet unknown coordinate which 
moves with the diffusive wave hj r = s{t) + z. The wavefront is then determined by s = 1 — z/, an equation 
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which is derived from the leading order terms of ( p.20| ) by taking expectation values. The speed of propagation 
is reduced slightly by the noisy coefficients, although this is a first-order effect, only being present in the 
(9(z/) terms. The large-time and large-size asymptotic solution is 



Crit) ~ Ic^erfc I ''^^ 1 , (3.21) 

' /2(l + z/)t; 

showing that the noise also broadens the diffusive wave in aggregation space as in Case I. The erfc shape is 
determined by the higher-order terms of (|3.20| ). 

The domination of aggregation over fragmentation in the contracted form of Case III is not altered by 
the presence of small amplitude noise 

An+l— 1 A„fl— 1 

An= E Sr + Oiu^), En+1= J] ^-+1- (3-22) 

r=A„ r=A„ 

In this case = 1 + A„ and the noise in the aggregation term, A„, is 0{h') whereas the noise in the 
fragmentation term, (3n+i = En+i, is much smaller, being of magnitude O^u^). Superficially this gives 
the impression of the contracted system being more strongly aggregation-dominated than the full system; 
however, this accentuated dominance is correct for similar reasons to those expounded for Case I. Since 
E[6r] = and Y[6r] = 0{v'^), we have E[A„] = and V[A„] = 0{v^Y also E[e,] = v and V[e,] = 0{v'^) 
implies E[K+i] = and V[K+i] = 0{y'^^). Equations (^[12]) and (g]^) imply 

^A^- i+EA"+o(^^) = n — 1+ ^ ^r^o[v^)\, (3.23) 

Vn=l ^"+1/ V n=l / \t=\ ^r+l) \ r=l j 

SO the contraction procedure does not lose information from the first two terms of the partition function-to 
see this, compare Tjv in the above with Qr in (|3.17|) with R = An- However, in this case the partition 



function does not play such an important role as in Cases I and II where it determines the large-time 
asymptotic solution which is approached; here, it is the steady-state solution which determines the large- 
time asymptotics. Since Er = 0{i>^), whilst A^ = O^u), perturbations to the forward coefficients infiuence 
the first correction term whereas those to the backward rates do not. The steady-state solution is thus 
modified to 

xT = 1 + Ai - A, + 0(z/2) = 1 + ^(5, _ S(^n-i)x+k) + 0{u^), (3.24) 

k=l 

which has constant fiux L = l+Ai + (9(z/^). The leading-order terms (x„ = 1, L = 1) agree with the full model 
and with the model without noisy coefficients, however the first correction term is not in agreement with the 
full solution of the noisy model. The contracted model predicts a steady-state fiux of L ~ 1 + Z]fc=i Sk + 0{u'^), 
whereas the full model has fiux J ~ 1 + (5i — £2 + C?(i^^), so agreement is limited to the leading order terms 
(J = L = 1) only, with the first-order correction terms differing, being Olu). In this case microscopic detail 
in the first correction is involved in determining the steady-state fiux, but the information is lost in the 
coarse-graining contraction, so the procedure only gives the correct result to leading order in u. 

To determine the kinetics of approach to steady-state in the reduced model we transform from x„(t) to 
^„(t) by Xn(t) = x^'^^ipnit) to gain 

1pn = (1 + A„)(^/'„+i - ^tpn) + Enilpn-l " V^n) + -^(^'n-l " ^n+l)- (3.25) 
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Since A„ = 0{i') and En = 0{i' ), the two-term continuum expansion of ( |3.25| ) correct to 0{i') is 

Thus, including the O(z^) terms, the diffusive wave has an expected speed of unity since E[A„] = 0. This 



differs shghtly from the result for the full model presented in (|3.21|) , where the nonzero 0{h') perturbations 
to the fragmentation rates caused the wave to be slowed. After taking expectations, ( p.26| ) is solved by 
ijj = |erfc((n — t)/\/2i) which yields the large-time solution 

Xnit) ~ i (1 + Ai - A„) erfc (^^^ ■ (3-27) 

Thus coarse-graining has altered the 0{v) correction terms in the kinetics of the approach to the steady-state, 
as well as the 0{v) corrections to the steady-state itself. However, the leading-order behaviour is faithfully 
reproduced. 

3.4 Case IV 

This generic class of behaviour describes systems dominated by fragmentation, in which larger clusters break 
up at a much faster rate than smaller ones. 

When noise is absent, the partition function is identically zero (for r > 2); however, the presence of noise 
alters the partition function to 

Thus Qr rapidly decays with increasing r since as r ^ oo, = (9(z/^'~^) and p > 0. The presence of 
perturbations to the rate coefficients provides a non-zero equilibrium solution Cr = Qr given by equation 
( |3.28D which is approached, as we shall now show, via a diffusive wave. Transforming to new variables 
iprit) = Cr(t)/Qr, we find 

= 6r{A+l - ^r) + [{r-lY + Er] {^r-l-^Pr), (3.29) 

where Sr > with E,[6r] = v for all r, and '¥\e,\ = for all r, with Sr taking both positive and negative 
values. In the large r and large time limit, ipr becomes smooth in r so it is valid to take the continuum limit 
which formally yields 

^ = Uir-ir + Sr + er]^-[{r-iy + er^5r]^. (3.30) 

This type of equation has a diffusive wave solution; we denote its position by s{t) and transform to a 
coordinate which moves with the wave, via r = s{t) + z. At leading order, this formally yields the equation 
s = {s — iy + Eg — Ss- Since E[£fc] = and K[6k] = u > 0, there is an C(i^) term in the equation for the 
expected value of s{t), namely 

s = (s-iy-u. (3.31) 

Thus the noise has a slowing influence on the wave, but since s ^ oo as t ^ oo, the leading-order term is 
So = Sg, and we have 

s,it) = [{l-p){t-to)Y/^'-P\ (3.32) 
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An 0{i') correction term can be calculated by putting s{t) = so{t) + z/si(t) into equation ( |3.31| ); one finds 
So(t) is determined by ( p.32|) and Si{t) by si = —1 +PSq~^si, which yields 

s, = ~^^7/^^':''^ + K{t - t,Yl^'-^\ (3.33) 
(1 - 2p) 

for some constant K. Thus, for p < 1, the wave experiences a deceleration which reduces the speed by a 
constant amount but which is insufficient to stop the wave since the growth in equation ( pj.32D has the faster 
growth rate in the limit t ^ oo. A higher-order effect is that noise broadens the wave very slightly, since noise 
in ( 3.3(J| ) has the effect of increasing the diffusion coefficient. The effect of this, however, is minimal since the 



diffusion constant which the wave experiences grows without bound in the large time limit. In cases where 
p < \, higher order terms from equation (|3.30| ) yield the shape of the front as = \eTic{z/ ^{2s{t) + 4i^t)) 
so that the large-time asymptotic solution (including the first correction term due to v) is given by 

Cr{t) ~ iQ.erfc ( /""^^^ I , (3.34) 
l2s{t)+Avtj 



where s{t) = So(t) + usiit) is given by equations (|3.32|) and ( |3.33| ). As t — +oo, the effects of perturbations 



to the rate coefficients decrease in significance, both in the position and in the shape of the wavefront. 

The rates in the contracted model are given by a„ = A„ and Pn+i = ''^^^ + En+i where the perturbations 
are determined by 

An+l— 1 A^l— 1 

An = A-^^ U ^r = Oii^'), = E ^ + ^(^')' (3-35) 

r=A„ r=A„ 



where in the latter, we have again made the approximation ( |1.13| ) valid for large r. Note that, as in Case I, 
ar = O^u^) whereas Pr+i = + C^i^): suggesting that while the contracted model equations have the same 
structure as the full ones, they are more heavily fragmentation dominated, due to the rate at which a cluster 
Xn+i = ca„+a grows from the cluster Xn = ca„ is much less than for Cr+i growing from c^. If we assume that 
the perturbations are independent and randomly distributed with 

E[Sr]=iy, E[Er+i]=0, Y[Sr]=Vsiy\ Y[er+i] = Vy, (3.36) 

then the perturbations to the contracted rates satisfy 

E[A„] = \-~P^u\ E[E„+i] = 0, 

V[A„] = \-^p^[{Vs + iy-iy\ V[K+i] ~ KAi-2Pn2p(^-i)z,2^ ^^■'^'> 

where the last formula includes a simplifying approximation valid for large r and large A, (namely replacing 
X]r=Al~^^^ with J^^^^~^ r^dr, to yield X^~^r~^). The formulae in ( |3.37| ) show that the contraction procedure 
makes the amplitude of the perturbations to the fragmentation rate dependent on the aggregation number, 
with larger sizes having a greater variance in the fragmentation rates. However, at larger sizes the standard 
deviation of the noise grows with 0{nP^^-^'^), thus the noise never becomes as large as the deterministic 
component of the rate (/3„ = (9(n^'^)). 

We now use this knowledge of the rates in the contracted system including the leading-order perturbations 
to construct a partition function for the contracted system T„ and, for consistency, verify that T„ = Q,. when 
n = Ar. Using ( pl2D and (|335| ), we find 

N-l N-l p 

log pAAr(log AT - 1) + ^ log A„ - E 



n=l n=l 



Ajv-1 A]v-1 

-pAAr(log(AAr) - 1) + ^ ^^^^r - E (3-38) 

r=l r=l ^ 
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where a similar approximation to that of ( |1.13| ) has been made. This result should be compared with Qp> for 
i? = Ajv from ( p.28|) , which yields 



R-l 



R-l 



logQn ~ -pR{\ogR- + ^og6r - E 



(3.39) 



Thus there are differences between and T^v, but at large these scale with logA^, which is of lower 
order of magnitude than the first two terms of T^r or Qajv- large N, the leading-order and first correction 
terms in the log of the partition function grow with iVlogA^ and N respectively, so both are reproduced 
correctly in the contracted description, as are the perturbations 6k,ek+i as can be seen by comparing (|3.38|) 
with ( |3.39D ; it is only high-order correction terms which differ. An alternative comparison can be made 
between 

((N-lV)p^ ((R-lVY 
mN]= ,}i,J-L-u E[Qh] = (3.40) 



X-pX{N-l),^X{N-l) - ~ - 

again showing that for large A^ and i? = A^v the leading order terms agree, since Stirling's formula leads to 
the dominant terms in both expressions being {\NY^^ e~^^^ / v^^ . 



In the present Case, the system tends to the modified equilibrium solution Xn = T„ given by ( p.38| ). The 
transformation ipnit) = Xn{t)/Tn enables us to find the large-time asymptotics; when applied to the equation 



Xri 



A72— 1X72— 1 I^n-^r, 



[n- 



Xji ~\~ 



n 



n+l 



it yields 



(3.41) 



(3.42) 



Since A„ = 0{v^) whilst En = 0{v) we keep only those terms involving perturbations to the fragmentation 
rate {En)] thus on formally taking the continuum limit, we find 



dt 2 



ly^ + En 



[{n 



(3.43) 



The substitution from independent variable n io z = n — s{t) yields the expression s = (s — l)^'^ since 
E[ii^„] = 0; hence the noise has no 0{i>) effect on the expected speed of the wave in the large-time limit unlike 
the full model (see equation ( p.31|) ). Thus as s — > 00, we have s = s^^^ and so s{t) = [{l—p\){t — to)Y^^^^''''^^ 
as in the full model with noiseless rate coefficients. As in Case I, the coarse-graining process has reduced the 
effect of the perturbations. These solutions are only valid for pX < 1; when pX > 1 the system instantly gels 
as noted by Brilliantov & Krapivsky and the large time asymptotics differ significantly. However, this 
case is less relevant physically, since the aggregation rate cannot normally grow faster than the cluster size, 
and usually grows much more slowly, giving rise to exponents strictly less than unity (i.e. p < 1). 



3.5 Case V 

This is perhaps the most interesting of all the nine classes of generic behaviour. Both aggregation and 
fragmentation are present and finely balanced; both occur faster at larger aggregation numbers than smaller 
ones. As expected in general crystal growth and dissolution processes, there is relatively slow nucleation of 
critical nuclei from the free "monomer" phase and faster growth/dissolution of supercritical clusters. 

In the case with no noise, aggregation and fragmentation are perfectly balanced, implying the same value 
for the partition function Qr = 1 for all r. The presence of perturbations to the rates alters this, to 

a = l + E(^^^^)+0(-')- (3-44) 
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Thus for z/ <^ 1 the leading order behaviour {Qr = 1) is unaltered, but suffers an O(i^) correction term. 
However, for p < 1 a result similar to that derived in Case II holds, where the effect of the perturbations 
accumulates, so that at large-r the noise may become a leading-order effect. For Case V we assume each 
6k, Ek has a mean of zero, a variance of z/^ and all are independent random variables. Thus at large r, 
has an expected value of unity with a variance of z/^r ^-P/{l-p). Thus when r = C»(z/-2/(i-p)), the variance 
is 0{1) and so the perturbations influence the leading-order term in the asymptotic expression for Qr. Note 
that if p > 1 then the variance of Qr approaches v'^Civ) (where C,{z) is the Riemann zeta function which 
arises in the solution of the noiseless time evolution for Case V; see ||TT|, ^), and so the accumulation does 
not become a leading order effect. 

In the case p > 1, however, the system does not evolve towards the equilibrium solution, but instead is 
attracted to a steady-state solution with a more rapid decay in the limit r — oo which we now examine in 
more detail. Perturbing the rate coefficients modifies this state from one of constant flux with J = l/([p) to 
J = l/C{p) + Ji, where 

1 ^ f 6k ^ 1 




" C(p)%?i i S 5: .J- (3-45) 



This gives the steady-state concentrations 




(3.46) 



Unlike the steady-states in Cases III and VI, the perturbations 61,82 do not play a special role in this solution; 
rather all perturbations influence the steady-state flux and concentrations. 

Let us now analyse the approach to equilibrium for p < 1. The equilibrium solution is identical to the 
partition function Qr since Ci = 1. We transform to new independent variables iprit) = Cr{t)/Qr, with Qr 
modified by the presence of noise, as in ( ^.44| ); for large r, we formally obtain the continuum approximation 



1 dip 
rP dt 




(3.47) 



This describes physical behaviour that approaches equilibrium neither by a diffusive wave nor by a simple 
diffusive process as Case II did; instead there is a more complicated similarity solution. Including O^u) 
correction terms, taking the expectation of this leads to ipt = {f^'ipr)r so the similarity solution 

_ .■--exp(-„^-7(2-ri^)d« 

fS°u-'exp{-u^-<'/(2-pr-)du • ^ ' 

provides the correct asymptotic approximation in the limit z^ <^ 1; the correction terms being much smaller 
than O(z^) (that is o(z/)), in contrast with Cases IV and VI where the correction terms are O(z^). 

When p > 1 a steady-state with non-zero flux is approached. We convert to new variables Cr{t) = cf^ipr{t)- 
For 1 < p < 2 the continuum approximation valid at large r is formally 

^ = (l I ^r + sA d'^jj ^ / 2 - p ^ 6r-er \dip 

rP dt \ 2rP J dr"^ \ r rP J dr' K ■ ) 

which, when we take the expected value of each term, reduces to 

dip pd'^ip , . p_idip 



+ (2-py-'^, (3.50) 



dt dr"^ dr 
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the correction terms due to 6r and Sr again being much smaller than O(z^), that is of magnitude o(z/). This 
equation has the solution 

/o~MP-2exp(-u2-p/(2-p)2)du ■ ^ ' 

Because p > 1, the perturbations are uniformly small in r in both advection and diffusion terms. Thus 
this is a uniformly- valid leading-order solution. For p > 2, the similarity solution of equation ( p.51 ) is not 
well-defined, and so does not determine convergence to the steady-state. 

In the contracted system of equations, the rates are given by a„ = n^^ + A„ and Pn+i = + En+i, 
where, again using (|1.13|) , we find 

^n = n^'j:t ^n.i = n^^E^- (3.52) 

r-=A„ r=A„ 

In contrast to Cases IV and VI, perturbations to the forward and backward rates in the contracted equations 
have the same order of magnitude. If the noise in the full model is assumed to be randomly distributed 
according to 

E[Sr] = 0, E[6r+i] = 0, Y[Sr] = Vsv\ V[£,+i] = V^v^ (3.53) 
then in the contracted model we have 

(3.54) 

As in equation ( p.371 ), the approximations replace I]r=An~^ by J^J^j^^~^ dr . From ( p.54| ) we see that the 
amplitude of the noise is size-dependent in the contracted description but, even at large aggregation numbers, 
it never rivals the deterministic component of the rate coefficients. The expected values and the orders of 
magnitude (in u) of the variance are correctly maintained in the reduced model. 

Following a similar method to that of Case IV, we insert the contracted rates (|3.52| ) into ( p.l2| ) to find 
the contracted partition function 

Ta. = 1+E ,a +^(^^) = 1+ E ^^-^ + 0(^^) (3.55) 

In this case the partition function constructed from the reduced rate coefficients matches the first two terms 
of the full (microscopic) partition function Qf{ with i? = A^v as in ( |3.44| ). This function also provides the 
large-time asymptotic approximation to the equilibrium solution for p\ <1. 

In cases where p > 1/A, a steady-state is approached instead of the equilibrium. This differs from the 
full model, where the transition from equilibrium to steady-state occurs at p = 1. Values of p satisfying 
1/A < p < 1 will approach a steady-state in the contracted description but the equilibrium solution in the 
full model. Thus when < p < 1, in order to capture the correct qualitative behaviour, one must ensure 
that A < 1/p. 

The steady-state solution for p > 1/A has flux L = l/({pX) + Li where 

the steady-state being given by 

" apX)tnk'' CiP^)rh\ rP^ ,i^_,,kP^ rP^t^kP' ' ^' ^ 
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This clearly has a similar form to (|3.46| ); however, whilst ( |3.57| ) displays the same qualitative behaviour as 



(|3.46|) , there are quantitative differences. The leading-order part of (|3.46| ) at large r decays according to 
~ 1 / {p — l)C,{p)rP~^ ^ implying that for r = A„ the quantity asymptotes to 1 / {p—l)C,{p)\P^^nP^^ whilst 
Xn ~ 1/ {p\ — l)C,{pX)nP'^~^\ while both decay algebraically, the coarse-graining procedure has altered the 
exponent of the decay. 

Let us turn now to the kinetics of the system's approach to equilibrium or a steady-state. Since both A„ 
and En are 0{i>) the kinetics of the contracted system is very similar to that of the full system. For p <l/\ 
the system approaches the equilibrium solution x„ = T„, where to 0{v), T„ is identical to Qj. with r = A„ 
- compare equations (|3.44|) and ( p.55|) . The manner of approach is remarkably similar to (|3.471 ), namely for 
ilj{n,t) = Xn{t)/Tn we formally have 

(3.58) 



which, on taking expectations leads to the solution 



I, _ /.7,v(.-..)^-^"exp(-^-^V(pA)^)rf^ 

' J,^uP^-^exp{-u-P^/{pXy)du ■ ^ ^ ^ 



When p > l/\ the system approaches the perturbed steady-state solution x^^^ given in ( p.56|) -( p.57|) according 
to 



dt 2nP^ )dn^^[ n ^ nP^ ) dn' ^^'^^^ 



where ip{n,t) = x„(t)/x^^^. Equation ( |3.6(J| ) corresponds to the kinetics of the full model ( ^.49| ). As in the 
full model, there are also corrections to the advection term which are more significant than the n~P^ 

terms if p > 2/A. After taking the expectation of equation (|3.60|) , the solution 

_ 47,i/(.-..) UP'-' exp(-^^-^V(2 - p\y) du 
^ /o°°MP^-2exp(-«2-pA/(2-pA)2)dM ' ^ ^ 

can be found. 



3.6 Case VI 

In this generic case there is virtually no fragmentation; the system is dominated by aggregation which occurs 
more rapidly for large cluster sizes. This scenario is typical of nucleation and growth processes in heavily 
supersaturated solutions, in which critical nuclei form relatively slowly, but then grow in size very rapidly. 

As in Case III, the partition function is not defined when noise is absent since all the fragmentation rates 
are then zero; when the fragmentation rates are perturbed by noise the partition function Qr can be defined 
but is sensitive to the amplitude of Sk 

thus Qr = O^iy-^^-^^). However in this case the system does not tend to the equilibrium state = Qr', 
instead, it tends to a steady-state given by constant flux J = 1 + (5i — 2~Pe2) + 0{h'''), which implies 



^sss _ J_ 



1+ 5^- 2-%2 ■ 



(r + l)P 



(3.63) 
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Thus, altering a,., 6^+1 by O(z^) alters the steady-state by an 0(z/) amount. As one might expect, the 
perturbations to rates at larger aggregation numbers are less significant than those to lower aggregation 
numbers; also as in Case III, the first two perturbations 61,62 influence all other concentrations at 0{v). 
Whilst the presence of noise restores the existence of the partition function, equation ( |3.63|) is the more 
physically relevant result, in that the partition function which increases with r will not be directly manifest in 
a simulation or experiment, where as the steady-state concentrations will be. These decrease with increasing 
r since p > 0. 

The manner in which the steady-state solution is approached is found by substituting Cr{t) = cf'^iprit) 
into the determining equations (|2.5| ) which, with 6 = and a = 1, gives 



Ipr = SriiJr-l " A) + (^^ + Sr){i'r+1 " Ipr) + "^(^r--! " A+l) , (3-64) 

where J is the steady-state flux into the system. We have two-term asymptotic expansions for J and cf^ given 
in, and just before, ( |3.63|) which enable us to simplify this equation. In the large-time and large aggregation 
number limits, we formally take the continuum limit 

% = + ^r + er)^~ (r- + Sr + er- 2..-,0 ^. (3.65) 



As noted earlier (Case IV, eqn ( |3.4| )), this equation has a solution in the form of an advective-diffusive wave; 
to determine the speed and shape of such a wave we transform to a frame of reference which moves with the 
wave by r = s{t) +z where r = s{t) denotes the position of the wave. The leading-order terms in ( |3.65| ) yield 



s = sP + 5s + es- 2es+i, (3.66) 

as s oo. Since E[5fe] = and 'E.[ek] = u > 0, taking expectations leads to s = ~ u, so perturbations have 
a small slowing effect on the progress of the diffusive wave. We solve ( |3.66| ) by assuming s(t) = So(t) + z/Si(t) 
with leading order solution ( p.32|) , as in Case IV. As Sq — >• oo, the 0{u) correction term is determined by 



Si = —1 + psq Si with solution ( |3.33| ). A higher-order effect is that noise broadens the wave very slightly, 
since noise in ( |3.65|) has the effect of increasing the diffusion coefficient. The effect of this, however, is minimal 



since the diffusion constant which the wave experiences grows without bound in the large time limit. When 
p < \ the shape of the wavefront can be calculated in the same way as for Case IV, the solution for ip{r,t) 
being identical, and leading to 

Crit) ~ ic^erfc I 1 , (3.67) 




where s{t) = So{t) + z/si(t) is given by eqns ( p.32| ) and 

In this case, coarse-graining maps the rates to a„ = n^^ + A„, f3n+i = En+i where the perturbations to 
the reaction rates A„, En+i are given by the formulae 

Artfl— 1 r Artfl— 1 

^n = nP^ ^ + 0(z/2), En+^ = A-P" n ^r+i = 0{u^). (3.68) 

r-=A„ ^ r=An 

Thus the domination of aggregation over fragmentation persists. We assume that the perturbations in the 
full microscopic model satisfy 

E[Sr] = 0, E[er+i]=u, Y[6r] = Vsu\ Y[er+i] = Vy, (3.69) 
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whence we find the perturbations in the reduced description satisfy 



E[A„] = 0, E[En+i] = A-pV\ V[A 



VsX 



l-2p^2p{A-l),^2 



A-2p^[(l/+l)^_l], 



.2A 



(3.70) 



These results are identical to Case IV (|3.37|) , but with aggregation and fragmentation rates reversed. 
Following the method used in the earlier cases, we insert ( p.68| ) into ( |2.12| ) and find 



Af-l 



n=l 

Ajv-1 



n=l 



logT^ ~ pXN {log N- 1) - E log^n+i + E 
~ pAiV(log(AA^) - 1 
which agrees with a direct expansion of Qr (|3.62|) 



k 



A]v-1 ^ 



(3.71) 



k=l 



k=l 



R-1 R-l r 

Ok. 



logQij ~ pR{\ogR- 1) - E log^fc+i + E 



k=l 



k=l 



kP' 



(3.72) 



For large R = An, the differences between log Qr and log T n grow with log R, or log N, whereas the quantities 
themselves grow at the much faster rates of RlogR and A^logA^. An alternative comparison can be made 
between 



N 



^-p\{N-l)^X{N-l) 

(iV- 1)!p^ 



and 



E[Qi 



,R-1 



{R-iy.p' 



(3.73) 



for R = Atv, where Stirling's approximation can be used to show that both expressions are dominated by 

However, in this case the system approaches a steady-state solution. When the leading order effects of 
the perturbations are incorporated, the steady-state becomes 



1 



1 + (Ai 



n 



(3.74) 



which has constant flux L = 1 + Ai. This has the correct qualitative behaviour, decaying as n increases. 
The leading-order term for large n is x„ ~ n~P'^ whilst the microscopic model predicts ~ r~P so that for 
maximum accuracy we should have Cr ~ {\n)~P^ when r = A„ which has a qualitatively similar shape (the 
decay is algebraic), but differs quantitatively since the exponent differs. The steady-state flux also differs. In 
the full model it is given by J = 1 + 5i — £2 which depends on noise in the fragmentation eate coefficient (£2) 
and has the expected value E[J] = 1 — z/ (since K[Sr] = and K[er] = v for all r). In the contracted model 
the flux depends only on the perturbations to the aggregation rates, with E[L] = 1 + o(v) since L = 1 + Ai 
and Ai involves the perturbations to the aggregation rates 5^ for r = 1, 2, . . . , A. 

To 0{y), the equations which determine the kinetics of the approach to steady-state are 



[n- 



'1)"" + A„_i 



X-n—l 



n 



+ A„, 



(3.75) 



since A„ = C(z/) and En = Oi^u^). The approach can be elucidated by substituting x„(t) = x^^^ipnit) into 
( |3.75D , where x^^** is as given in p.74| ), since then 



= {nP^ + A„)(^,+i - + —ii^- 



n— 1 



n+l; 



(3.76) 
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The steady-state solution x^^ has as an asymptotic approximation (p.74|) with steady flux L = 1 + Ai. In 



the large-n and large-time limit, we take the continuum limit, obtaining the partial differential equation 

Following earlier analysis, we transform to a moving coordinate frame by n = s{t) + z, obtaining the equation 
s = s^'*' + Ag for the position of the wavefront. Since ]E[As] = 0, there is no effect of the noise on the 
leading-order progression of the wave through the system, and we have 

s{t) = [{l-p\){t-to)Y/^'-^'\ (3.78) 

the corrections due to noise being of order O^u"^) and above. For pA < | the shape of the wave is given by 
il) = |erfc(zA/l — 2pX/ J2s{t)), so that the large-time asymptotics are governed by 



2s{t) 

for n — s(t) ~ \Js{t) as t ^ oo. For | < pA < 1 the shape of the wave depends on the initial data for all 
time. These results are consistent with Cases I, III, and IV, where the coarse-graining procedure reduces the 
effect of perturbations to fragmentation rates. 

3.7 Case VII 

Here fragmentation is the dominant process, with smaller clusters shedding monomers much more readily 
than larger clusters. 

The coagulation and fragmentation rates are given by = Sr and fe^+i = + ^r+i with p < 0, and 
in all of Cases VII, VIII and IX, the perturbations are taken to be positive random parameters distributed 
according to 

E[6r] = z/, E[er+i] = ly, Y[6r] = Vsi^^ V[e,+i] = V^. (3.80) 

As with Cases I and IV, if there are no perturbations to the rates then the partition function is identically 
zero; when noise is introduced we find 

«-n^. (3.81) 

so that Qr decreases rapidly in magnitude as r increases {Qr = 0{i>^~^)). However, since p < 0, this ceases to 
be valid when r = (9(z/^/p), because the perturbations have the same order of magnitude as the deterministic 
part of the coefficients. When rz/^/^^~^) < 1 the approximation (|3.28| ) is valid; however, this ceases to be hold 
when r = (9(z/^^/*^^~^)), where the variance of the sum reaches 0{l). At this point the random perturbations 
affect the leading order behaviour of the partition function. 

At large times the system approaches its equilibrium solution. As in Case IV, the equation determining 
the progression of the diffusive wave is formally s = {s — iy + Es — Ss- On taking expectations of this equation, 
all i/-dependence disappears, indicating that the noise has no net effect on the progress of the wave. The 
equation s = is valid until s = 0{u^/^)^ hence s{t) follows ( p.32| ) until t = 0{u"^^^/^). After this time, the 



random components of the coefficients influence the leading-order motion of the wave. Thus there are subtle 
differences between Case IV (p > 0), where ( p.32|) holds for all time, and Case VII {p < 0), where (|3.32| ) holds 



only for times up to (9(z/^/^). At r = 0{v^/^) the random component of the rate coefficients in the diffusive 
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term in ( |3.30| ) becomes leading-order. Since this term contains the sum of two positive perturbations, there 



is a net increase in the diffusivity, hence the wavefront widens as described by ( p.34|) . 



Following the coarse-graining contraction, formulae ( p. 35 ) still hold, so the noise satisfies 



A„+i — 1 -j^ 



E[A„] = A-PV^ E[K+i] = n^V Yp 



k=Ar, 

An+1-1 



(3.82) 



fe=A„ 



Due to the decay in fragmentation rates at large r, there is an aggregation size at which the perturbations 
assume the same order of magnitude as the deterministic part of the rate coefficients. To find this threshold, 
we equate n''''^ with the expectation of -En+i in ( p.82|) . This yields n = z/^/^A^^^^/^, corresponding to an 
aggregation number of r = u^^^X^^p in the full description of the model (before contraction), showing the 
correct order of magnitude when compared with the full model, which gives r = z/^/^. 

Using the coarse-grained rate coefficients to construct a partition function Tat, we find ( |3.38D holds for 
N = C(l). Together with equation dS^S^) and ([OHD , for R, N = 0(1), these imply 

^l^"l=(iVrT)Sx(l-^i:;^) ElQ«l = ^^(l-.|:-]. (3.83) 

so that for i? = A^v we have both expressions growing with R (or equivalently with A^), whilst their difference 
only grows with log R. This shows nice agreement . At large R, N the formulae ( p.83| ) fail due to the noise in 
the fragmentation rate becoming leading order. This occurs for R = 0(z/^/(P-^)) in the latter case and, upon 
usmg ( p:82D , A^ = C(z/i/(P-i)) in the former; again showing good agreement between the full and contracted 



models. 

As with Case IV, the system converges to equilibrium via a diffusive wave which travels from small to 
large-n, invading the region where the initial conditions x„ = have not yet been modified, and leaving 
behind (at smaller n) the equilibrium solution x„ = T„. Thus we use the substitution ipn(t) = a;„(t)/T„, 
which satisfies {^^. Since E[A„] = 0{u^) and E[E„+i] ~ A^-Pn^^^-^)// from j ^M l, the continuum limit of 



(|]42D to 0(z/) reduces to 



At large-times, the leading-order term determines the speed of the wavefront (s), giving the equation s = 
gP^ _l_ X^-PgP'^'^^^^i'. As z/ 0, the solution of this equation can be formulated as s{t) = so{t) + z/si(t) where 

so(t)= [(l-pA)t]^/(i-^^), and ,(t) ^ ,o(t) + . (3.85) 

This shows that as the wave reaches larger values of n, perturbations to the rate coefficients infiuence the 
leading- ordermotion of the wave. This occurs when s ~ n = 0{h'^^^^~^^), or when t = 0{u~^~^^^~^^^^p~^^), 
which is the same order of magnitude as in the uncontracted case. At this order of magnitude the perturba- 
tions to the rate coefficients also cause the shape of the wavefront to be modified. 

3.8 Case VIII 

Here aggregation and fragmentation are finely balanced, but occur much more rapidly at small cluster sizes 
than at large cluster sizes, a situation that may arise in growth processes in solution below saturation without 
nucleation barriers. 
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In the noise- free case this system converges to the equihbrium solution Cr = 1. When noisy coefficients 
are introduced, defined by = af" + 6r, K+i = br^ + e^+i, with p < and the perturbations distributed 
according to (|3.80|) , then the partition function is modified to 

and so for r = 0{1), the equihbrium solution can be approximated by 

Cr = l + E k-''{Sk - Sk+i) + 0{u^). (3.87) 

k=l 

However, this approximation ceases to be valid at large r due to the perturbations becoming a leading-order 
effect. Given that the variance of each perturbation is (9(z/^), the variance of Cr in equation ( 3.87|) is (9(z/^r^~^) 



so that at r = u~'^^^^~p\ the leading order expression for Cr is no longer unity. At this aggregation number 
the cumulative effect of the perturbations becomes as important as the leading-order term. 

When we turn to the analysis of the kinetics of the approach to equilibrium, we note the similarities 
between this case (where p < 0), and the approach to equilibrium in Case V for < p < 1. The continuum 
equation ( |3.47| ) is also valid here, but whereas in Case V E,[6r + Sr] = 0, here K[6r + Er] = 2z/, so when we 



take the expectation of (|3.47|) we obtain an equation with an increased diffusivity, namely 

Thus there is a large-r region where the perturbations influence the leading-order kinetics. This occurs 
for r = (9(i/^/^) and larger, where noise in the diffusion term is comparable with the deterministic part 
of the diffusivity. In the advection term, the noise has zero mean. Thus for rz/^^/^ < 1, we expect the 
similarity solution ip = f{r]) with i] = r/t^/*^^"^-' to be valid and il>{r^t) determined by equation ( p.48|) ; but 
for r > (9(z/^/^), the leading-order kinetics are influenced by the noise. 

In the coarse-grained contracted system, the forward and backward rate coefficients are given by 

An+l— 1 Artfl— 1 

«n= X{{rP + 5r), Pn+1= Ui^' + ^r+l), (3.89) 

r=A„ r=k„ 



and whilst ( p.52| ) is a valid approximation for small n, it breaks down at larger n. Assuming that the noise 



is distributed according to ( p.80|) before the contraction procedure, we flnd 



-^n + l — 1 1 A„-|-l — 1 



(3.90) 

following the coarse-grained renormalisation. The even balance of aggregation and fragmentation is reflected 
in the identical expectations of A„ and -E„+i and in the same order of magnitude of their variances. When 
n = Oi^u^^P) the formulaic component of the rates (O^u^'^)) matches the random element (providing A is not 
large); this aggregation number agrees with the analysis of the fully microscopic model. 

We now use the coarse-grained rates to calculate the partition function (T„) for the coarse-grained system. 
Equation ( |3.55| ) is valid for r = 0{1), however, since p is negative the sums are divergent at large n and large 
r. Although E[Tjv] = 1 = E[(5_r] for all N and for all R, the variance of T^v grows with N according to 

(Af-1)A+1 , 

Y[T^] = u\Vs + V,) E (3-91) 

n=l " 
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which agrees with V[(5_r] for R = A„. The asymptotic expression ( |3.55| ) thus ceases to be uniformly vahd 
when the variance becomes 0{1), namely when n = giving the same order of magnitude of 

aggregation number (r = A„) as can be derived from the full system of equations (|3.44| ). 

The kinetics of approach to equilibrium are similar to Case V with p < 1/A; that is by a perturbed 
similarity solution rather than a moving diffusive wavefront. Since in the present case the perturbations obey 



equation (|3.80|), when we take the expectation of equation (|3.58|) we obtain 



„-.^^ = (l + 2A->',rV)^| + P^^. (3.92) 

at ^ ' aw^ n on 

Thus when one ignores the 0{i') terms, a similarity solution can be found; and for n < 0{u^^^), we expect 
the similarity solution to give the leading order behaviour. This solution has the form tpn{t) = f{v) with 
r] = n / / as in ( p.59[ ). However, for n = 0{u^^^) and larger, another more complicated solution takes 



effect, where the noise influences the dynamical behaviour at leading order. For n larger than 0{u^^^) the 
system behaves as if all rates had been chosen at random. 



3.9 Case IX 



The final case corresponds to a system in which aggregation dominates fragmentation at smaller cluster sizes; 
at larger cluster sizes, the perturbations to both aggregation and fragmentation rate coefficients entirely 
swamp this effect. 

Here, we have = + 6r, &r+i = £r+i with the perturbations Sr,Sr+i distributed according to (|3.80|). 



Formally the presence of noise allows the partition function to be written as 



Qr 



k=l ^k+l 



(3.93) 



With noiseless rate coefficients, the system evolves to the steady-state with unit flux ( J = 1); however, using 
the perturbed coefficients in equation (|1.5|), a calculation of the steady-state flux leads to J = since at 
large values of r, the term 1/arQr has a positive expectation value which is independent of r, making the 
sum divergent. Thus in this case the presence of noisy coefficients changes the large-time behaviour from an 
approach to steady-state with unit flux to an approach to the equilibrium solution, which is given by = Qr 
where, for r = Qr is given by equation (p.62|) . This expression, however, ceases to be valid for values 

of r larger than (9(z/~^/'^^^p)). At intermediate times, before mass has been transported to larger aggregation 
numbers, we expect the steady-state solution to be manifest. At larger times, the system undergoes a complex 
transition from a steady-state solution to the equilibrium solution. Evolution will then be dominated by the 
random perturbations to the rate coefficients, and will follow the kinetics described in the next section. 

In the coarse-grained description of the system the equation for A„ in ( |3.68D is valid for small n, but 
ceases to be valid when n = (9(z/^/^), where the noise becomes leading order; thus there is an intermediate 
asymptotic regime where the system remains aggregation-dominated. Using ( |3.8CI| ) we find the distribution 
of perturbations to be as follows: 



fc=A„ 

A„+i-l 
fc=A„ 



(3.94) 



however, at larger aggregation numbers the system is dominated by the random perturbations. 
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Calculating the partition function T^v from the coarse-grained rates we find ( p.71D holds for small and 



intermediate values of N. For large N however it fails since p is negative and the second term becomes the 
same order of magnitude as the first when N = 0{u^/p). The presence of noise transforms this case from one 
which approached steady-state solution (with L ^ 1) to one which approaches equilibrium (i.e. L = 0); as in 
the microscopic system. Thus the coarse-graining procedure has faithfully retained the structural difference 
that noise has made to the full system of equations. 

In the large-time limit, this case shows convergence to the equilibrium solution, since at large aggregation 
numbers, the system appears identical to the case with all rates chosen at random described below, in Section 
3.10| . In this case the kinetics are dominated by the random perturbations to the rates, so it is impossible to 



give a detailed description of the time-dependent solution. 

At intermediate times, when there is little mass in large aggregation numbers, however, we expect the 
system to behave as in Case VI, namely the approach to a steady-state solution. Since the rate coefficients 
at small and intermediate aggregation numbers are dominated by aggregation, we expect a diffusive wave to 
move into the larger-n region leaving behind the steady-state solution Xn = n~^^. Since this moves at a rate 
given by ( pTHD , when t = 0{u^/p-^) the wave reaches aggregation sizes n = (9(z/^/p), where the perturbations 



to the rates are of the same order of magnitude as the deterministic component of the rates. Thus after this 
time the system undergoes a transition from the state 



^-^'^ ~ forn » u^P (^-9^) 



to the equilibrium state x„ = 



3.10 The system with totally random rates 



The last three subsections show that, at large times and with the rate perturbations as defined in Section |2Tl| , 
in each of the Cases VII, VIII and IX the temporal evolution tends to the equilibrium solution rather than a 
steady-state solution; at large r the partition function Qr is dominated by the perturbations, hence so too is 
the equilibrium solution in all three cases. From this perspective, these three cases can then be thought of as 
lying in the same universality class, namely a Becker-Doring system with "totally random" rate coefficients. 
By this we mean a formulation of the Becker-Doring equations (|1.4|) in which all the rate coefficients a^, K 
are chosen at random (but remain independent of time); in the notation of this paper, we have a = b = 0, 
leaving ar = Sr, br = Er- Here the rates Sr, Sr are independent random parameters whose expectation and 
variance is size-independent. 

If we take the monomer concentration to be fixed at unity, the governing equations are 

= 5y—\Cy—\ — E^Cf — Sj.Cj. -\- £j._|_iCj.^i , r ^ 2. (3.96) 

It is possible to define a partition function from the coefficients 6r, Sr provided none are zero. If we assume 
that in the large-time limit the system will tend to a steady-state 



V k=l J 



then imposing the fastest possible decay on the concentrations in the limit of large r, we find 

1 oo 1 

7 = E TTT- (3-98) 
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Since neither Qk ^ oo or 6k oo as k ^ oo the terms in the sum do not decay, so the sum diverges. We 
thus have J = 0, and the system must tend to the equihbrium solution Cr = Qr- To find the kinetics in the 
large-time limit we substitute = Qri^r] then ipr satisifes the difference equation 

Ipr = — ^r'lpr — ^ri^r + '^rV'r+l, (3.99) 

This is a system of equations similar to that with which we started out in equation ( |3.96| ); the essential 
difference is that — >■ 1 as t — > cxo. We expect iprit) to be slowly varying in the large-time limit, in turn 
enabling us to take the continuum limit which gives the following partial differential equation for ip^-'- 

|.i(,,,,„,|| + „,_,,)^. (3,00) 

Thus the advective component of the process {5r — ^r)^ has zero mean (since ¥\5r] = ]E[e,.] = z/), whilst 
the diffusive component is always positive, having mean u. The expectation of the solution is thus ip = 
^eTk{r /2^/ui) and so we see that the system shares several similarities with Case II, notably convergence to 
equilibrium via the same kinetics albeit on a slower timescale. 

In Cases VII and IX we expect the waves to cease moving when they reach the large-r region where 
perturbations become dominant, namely when the wave front reaches r = (9(z/^/p); which occurs after a 
time of t = 1/(1 — . However, in all three Cases VII-IX the final approach to equilibrium is by a 

predominantly diffusive mechanism. 



4 Discussion 

A previous paper described in detail nine generic classes of behaviour into which the asymptotic dynamics 



of the Becker-Doring equations with power-law coefficients falls [T^. The nine classes capture qualitatively 
different physical properties which are shared by all models within the same class. In the present paper, we 
have concentrated on the detailed analysis of these nine cases; for each case we have considered the effect 
of the perturbations to the rate constants on the equilibrium or steady-state solution and on the large-time 
asymptotics. 

In Cases IV- VI, where the rate constants grow with increasing cluster size, larger perturbations could 
be considered, (that is the perturbations to the forward and backward rate coefficients 6r,er are each 0{1)). 
Because the unperturbed rates grow with r, there would then be a large-r region where the perturbations 
are insignificant and a smaller-r region where they should be taken into account. Since we are dealing 
with a discrete system, the small-r region consists of a complicated finite dimensional system joined to a 
simpler large-r system whose behaviour we already know by asymptotic analysis (King & Wattis 13). By 
taking a coarse-graining contraction with A sufficiently large, the whole of this small r regime could be 
mapped to a small dimensional system with only a few concentrations needing to be retained. As in the 
unperturbed problem analysed previously |TT], coarse-graining retains the correct leading-order structure 
of the problem, although critical exponents of p = 1 where the behaviour changes are mapped to p = 1/A. 

Cases VII-IX are much more complex, since here the rate coefficients decay with increasing aggregation 
number. Thus in the limit r — > oo, the perturbations 6r,er dominate the partition function (Qr) and the 
large-time kinetics. If we assume the perturbations are small (as we have done here), there is a small-r region 
where regular behaviour occurs, and a larger-r region where perturbations dominate the rate coefficients, 
the partition function and the kinetics; in this region the system behaves as if the rate coefficients had been 
chosen at random. Following an application of the coarse-graining process the positions of these regions 
remain invariant and, in the small r region, the leading order behaviour agrees with that of the full model. 
The only case in which the large time solution which is approached is altered by the presence of noise is Case 
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IX, which changes from convergence to a steady-state solution in the large-time limit into a system which 
converges to equilibrium. The equilibrium solution to which Cases VII and VIII approach in the large-time 
limit also suffer major modification due to the presence of noise. 

In renormalisation-theoretic jargon, Cases I- VI can be referred to as universality classes, whereas this 
is not in general true for Cases VII-IX, because the perturbations destabilise the steady-state solution, 
modify the equilibrium solution and alter the associated large-time dynamics of convergence to equilibrium. 
However all of Cases VII-IX share the same large aggregation number and large time behaviour, and so fall 
into a seventh Universality class, which we have referred to as the system with totally random rates (Section 
3.10| ). This classification occurs with the rate perturbations as described in Section |2.1| , namely those with 



size- independent mean and variance. 

If an alternative form of rate-perturbations were applied to the system, then an alternative universality 
classification would be produced. For example, if the rate-perturbations had the same size-dependence as the 
power-law component, that is E[5r] = i^r^ and E[£r+i] = i^r^ then the large-time asymptotics of Cases VII-IX 
suffer no leading-order modification due to the perturbations and each of these cases then corresponds to its 
own universality class. However, this one-to-one correspondance of our nine special cases with universality 
classes only holds in situations where the added noise decays faster than the specified power-law for the 
rate-coefficients as r —>■ oo. We believe such scenarios to be somewhat artificial, thus in this paper we 
have concentrated our analysis on the more generic and interesting case where the behaviour of the rate- 
perturbations at large aggregation numbers differs from that of the power-law. 

In future work, we plan to extend and generalise these renormalisation-theoretic results to a range of other 
Becker-Doring systems, including ones with different forms of rate coefficients which describe quite distinct 



physicochemical processes ^Uj , as well as the constant-density formulation of the Becker-Doring equations. It 
turns out that, for certain well motivated choices of rate coefficients, the renormalisation procedure is exact 
(in the sense that the approximation made after equation ( |1.13| ) is not necessary). The constant density 
Becker-Doring system has the added complication of being inherently nonlinear, and our renormalisation 
scheme accentuates the nonlinearity, making analysis of this problem even more challenging. 
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